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Abstract 

Grover's quantum algorithm improves any classical search algorithm. We 
show how random Gaussian noise at each step of the algorithm can be mod- 
elled easily because of the exact recursion formulas available for computing 
the quantum amplitude in Grover's algorithm. We study the algorithm's in- 
trinsic robustess when no quantum correction codes are used, and evaluate 
how much noise the algorithm can bear with, in terms of the size of the phone 
book and a desired probability of finding the correct result. The algorithm 
loses efficiency when noise is added, but does not slow down. We also study 
the maximal noise under which the iterated quantum algorithm is just as slow 
as the classical algorithm. In all cases, the width of the allowed noise scales 
with the size of the phone book as iV -2 / 3 . 
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I. INTRODUCTION 



There exist problems where the algorithm that solves them scales exponentially as the 
size of the input is increased, for example computing all possible chess games, factoring a very 
large number, etc. This dependence on the size makes them physically unsolvable for large 
enough inputs. Quantum algorithms have been invented to bypass this problem, like Shor's 
ffl] that turns tractable the problem of factoring numbers, and Grover's M that improves 
the classical search for an item in a phone book. In fact, the classical search algorithm does 
not scale exponentially. Rather, it is linear in the size of the phone book; Grover's quantum 
algorithm improves it to a square-root dependence. Recently, an experimental application 
of a quantum algorithm has been implemented fTj, and agreement between theory and 
experiment was found. 

Nevertheless, the strength of a quantum algorithm is also its weakness: a quantum 
computer performs simultaneous operations over large superpositions of states, which are 
very sensitive to decoherence. Fortunately, quantum correction codes have been developed 
with which a quantum computer can recover from errors in the presence of moderate 
decoherence. But these quantum correction codes are subject themselves to decoherence, 
and it is not fully understood how decoherence affects the correction itself. In this work, 
we study the intrinsic robustness of Grover's algorithm, when quantum correction codes are 
not implemented. 

II. GROVER'S QUANTUM SEARCH ALGORITHM 

Any classical algorithm for finding an item in a randomly ordered phone book (whether 
deterministic or probabilistic) requires N/2 steps on the average, because the only way to 
perform the search is to analyze each item one by one until the searched-for item is found. 
Recently, Grover invented a quantum algorithm that runs like 0(y/N). Let us review 
it briefly. 
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In a phone book with N = 2 n entries, each item can be represented by a binary label 
of length n or, equivalently, by a pure state of n spin 1/2 particles. The algorithm is based 
on constructing a coherent superposition of all these states, and applying repeatedly certain 
unitary transformations to it. 

Assume, for concreteness, that the item we are looking for is represented by the state 
III • • • |), i.e. by n spin-down particles. 

The algorithm works via the repeated action of the unitary steps below, starting from an 
initial state which we take to be the full coherent superposition of all states in the system, 
namely 



(2.1) 



V 1 / 

Of course, one could start equally well with some other initial state ||. 
The two unitary steps to be repeated are the following: 

First, invert the phase of the looked-for state trough the unitary transformation 



/ 



U-i 



-1 ■■■ 

1 : 

: '•• 

••• 1 



(2.2) 



Secondly, invert, with respect to the average, the phase of the looked-for state trough 
the unitary diffusion matrix 

2 
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These two steps are equivalent to the action of the following single unitary transformation: 
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When the unitary transformation U has been applied m times to the initial state \l/o, 
the new quantum state will be 



m m = U m ^ 
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The action of U on the initial state \&o yields only two distinct amplitudes A m and B m , 
whereby it is possible to recast the recursion relation in just two dimensions. The restriction 
of U to this two-dimensional subspace will be denoted by S. Explicitly, the amplitudes A m 
and B m are given by the recursion formula 
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The two-dimensional matrix S has eigenvalues e ±lip , with cosy? = 1 — 4, whereby 

1 
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thus 
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(2.8) 



'N \ ' VN^T 
From (|2.7|) , the probability of finding the state we are looking for if we measure ^ m is 



P(m) = \ A m \ = — [cos imp + \/ N — 1 sin my? J 



(2.9) 



With the change of variables y? = 26, P(m) can be written as 



P(m) = sin 2 (6 (2m + 1)) , (2.10) 
Clearly, P(m) is periodic, with maxima at 

9 (2m + 1) = mr, n integer, (2-11) 



The first maximum for large N is approximately at 

ttx/N 
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(2.12) 



and P max = -P(^imax) — 1- The number of steps required to find the state with almost 
certainty scales like yN, as shown in Q2.12 ). 



III. MODELLING NOISE IN GROVER'S ALGORITHM 

As stated in the introduction, quantum correction codes have been developed and it 
is supposed that in the presence of low but physically realistic levels of noise they are 
useful |5|,|J . These codes can be implemented only if a small enough subset of the quantum 
computer's q-bits undergo errors, and when the probability of occurrence of an error in the 
computation is lower than a certain bound. On the other hand, the real effect of the noise 
introduced by these correction codes over the original algorithm is not completely known, 
because they are quantum computations too. Hopefully, such errors are small and tractable. 
But what happens if they are not? Or, even worse, what happens if many q-bits undergo 
errors? Is it still possible to make sense of the computation under this hypothetical noisy 
situation when quantum correction codes do not suffice or cannot be implemented? If it 
does, how much noise the algorithm can bear with on its own? We now turn to the answer 
to these questions. 

In the particular case of Grover's algorithm, there is a simple way to model noise, because 



of the explicit recursion formula ( |2.6| ) for the amplitudes of the searched-for state. 

Suppose that in each step of the algorithm, a white or Gaussian noise modifies the state 
of the whole phone book according to 
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(3.1) 



where S is denned in (|2.6|) , and both a m and b m are noise, determined randomly by the 
standard deviation a (common to both, for simplicity) of their Gaussian distribution. Of 
course, the new state \P TO +i is appropriately normalized (that's what the denominator Norm 
is for). Explicitly, 



(a ^ 

' 'in 



-la log^i 



V 



(3.2) 



cos 2nx2 

where x\ and x<i are computer-generated random variables uniformly distributed over the 
interval [0, 1]. The two Gaussian variables a m and b m are mutually independent, and change, 
randomly, from one iteration of equation (|3.1|) to the next. Note that when a = 0, a m and 
b m are always zero and thus there is no noise. 

A crucial caveat is in order here: note that we introduce only two different errors, one for 
the searched-for state and one for all the other pure states. This approximation is physically 
unrealistic, but worthy of study. The full noisy situation would call for allowing iV different 
random variables to be added independently to each of the iV components of the state vector, 
instead of restricting ourselves to noise in the two-dimensional subspace where 5* (instead 
of U) acts. 

Now, we want to find the maximal allowed noise, quantified by a, in terms of both (a) the 
size iV of the phone book and (b) a given probability P cut for finding the searched-for state 
after a suitable number of iterations. If we set P cut = P m ax, then of course o can only be 
zero. As we allow for a decreased certainty of finding the result, and thus decrease P cut , the 
algorithm can bear with an increasing amount of noise. In the absurd limit of being happy 
with P cut ~ 0, which means we will not find the result, then any amount of noise is allowed. 
Of course, for any given P cut > 0, a large enough noise will destroy the algorithm. In the 
next section we establish the dependence of this maximal allowed noise, <7 max , in terms of 
and P cnt . 
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A. Computations and Results 



To find when the algorithm breaks down as we increase the noise, we treat the noise as a 
perturbation on the exact algorithm (recovered when a = 0). Beforehand, we fix the phone 
book's size N and the desired probability of finding the result, P cu t- 

First, we take a very small initial value of a and evolve the initial state \l/o i n equation 
fl2-lp according to the noisy iteration given in equation (|3.1| ). After m iterations, the prob- 
ability P(m) of finding the result is still | A m | 2 , where now the amplitude A m includes m 
additions of noise. It turns out that, on the average, P{m) still reaches its maximum after 
m max steps. This is a pleasant surprise. At first thought, one could have imagined that 
noise not only decreased P max (as it does), but also slowed down the algorithm (which it 
does not). To maximize the likelihood of finding the result we must measure the quantum 
state after m max iterations, with m max given by the noiseless equation ( p.!2| ). 

Now we compute P max = P(«Vax) and compare it with P cut . If P max is greater than P cut , 
we increase the value of a and repeat the computation, otherwise we stop (see the Appendix 
for details). In this way, we find the maximal a, labelled cr max , which is the limiting noise 
for P max > P cu t- Because of the probabilistic nature of the computations, we repeat this 
computation of a max many times (two-hundred): the value of a max we exhibit is the average, 
with a statistical error. 

We have carried out the evaluation of <r max for seven different phone book sizes N = 2 n 
(with n from 10 to 16) and for five different values of P cut (from 0.9 to 0.5 in steps of 0.1). 

For fixed P cut , the dependence of <T max on N is always of the form: 

ow(iV,P cut ) =a(P cut )Nt (3.3) 
where is a true constant, found to be 

= -0.696 ± 0.027 (3.4) 
and a varies smoothly from .9 to .15 as P cut decreases from .9 to .5 (see Appendix). 
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One of our main results is that the amount of noise that the algorithm can handle 
decreases roughly as iV~ 2 / 3 with the size N of the list. In general, since the number of steps 
needed in each iteration is of the order of N 1 ^ 2 , and at each step we add a noise of width a, 
we expect the maximal allowed cr max to decrease with N faster than N~ 1 ^ 2 . Equivalently, 
we expect <fi to be smaller than minus one-half. The actual value found, equation ( |3.4| ), 
satisfies this bound. We have not found a general analytic argument to pin down the actual 
value of 0. 

Alternatively, keeping A fixed instead of P cu t, the relation between cr max and P cnt can be 
written as 

ow (A, P cut ) = 7 (iV) - S (N) P cut (3.5) 

where 7 goes from 0.0024 to 0.00015, and 5 from -0.0020 to 0.00013 (log 2 iV = n = 10 and 
16, respectively), with errors of about 10% (see Appendix for details). This means that the 
width of the maximal white noise that may be allowed increases linearly with decreasing 

-Pent- 



Note that equations (|3.3| ) and Q3.5|) are just convenient slices of a surface in the three- 



dimensional space with co-ordinates (A, P cut , a n 



IV. GROVER'S ALGORITHM IS USEFUL EVEN IF P CUT < 0.5. 

In the derivation of the above results we exploited the experimental fact that the number 
of steps needed to find the searched-for state does not change when noise is present. Thus, 
another way to estimate the real maximal noise that the noisy Grover's algorithm can handle, 
while still improving the results of the classical search algorithm, is to let P cu t be even lower 
than 0.5. We now explain this. 

Since m max = jVN is always bigger than A/2, there is an integer In such that i7vm max < 
N/2, namely 

I N ~ -VN (4.1) 

7T 
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Therefore, we can repeat the quantum search In times with a low P cu t such that 1 — 
(1 — P C ut) lN > 0.5. We are assured that we will find the searched-for state with probability 
one-half in the same number of steps as the classical algorithm. Of course, the classical 
algorithm find the result for sure, and compared with that finding the result only half the 
time is not very satisfactory. Instead of .5, we could equally well have chosen some other 
(higher) probability to be satisfied with, but we take .5 for definiteness as the extreme, 
illustrative case. The point is that the P cn t we need to enforce on the noisy quantum 
algorithm is smaller than .5. Note also that we are disregarding the log 2 iV steps needed in 
each of the Jjv independent iterations to prepare the initial state ^o- Including them would 
of course lower a bit the maximal allowed noise. 

The limiting probability at maximum with which the iterated quantum algorithm is as 
slow as the classical one is 

Pcut>l-0.5^ 2 ^, (4.2) 

The meaning of this is, again, that we can let P cut be smaller than 0.5 for a given iV because 
if we run 7jy times the quantum algorithm with m max ~ iry/N/A steps, we will find the 
searched-for state with a probability of at least 0.5, and the total number of steps will be 
less or equal to N/2 (ignoring the log 2 N steps required for constructing the initial state 
*o)- 

To estimate this maximal noise that the quantum algorithm can bear before it slows 
down all the way to equivalence with the classical one, we proceed as follows. First, we 
choose the size N of the list to be searched, and keep it fixed. Then, using the bound 
(f4.2p, we determine P cut , which is very low. Finally, equation Q3.5D yields cr max , which is now 
significantly higher. For a variety of N, our results are shown in Table 1. 
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TABLES 



N 


Pcut 


"max 




1024 


0.034 


2.33 x 10~ 3 


1.0 x 10~ 4 


2048 


0.024 


1.48 x 10" 4 


4.3 x 10~ 5 


4096 


0.017 


9.03 x 10~ 4 


2.6 x 10" 5 


8192 


0.012 


5.68 x 10- 4 


1.6 x 10~ 5 


16384 


0.0085 


3.28 x 10~ 4 


1.7 x 10" 5 


32768 


0.0060 


2.13 x 10~ 4 


1.7 x 10" 5 


65536 


0.0043 


1.17 x 10~ 4 


1.1 x 10~ 5 



TABLE I. In the iterated quantum algorithm, for various sizes N of a phone book, the absolute 
maximal allowed Gaussian width <r max of the white noise, and its statistical uncertainty (between 
5 and 10 %). Also shown is the (low!) limiting probability P cut at maximum. 

In figure 1, we plot <r max as a function of N for the data of Table 1; the equation which 
fits it is 



a max = (0.275 ± 0.031) iV(-°- 68±001 ), 



(4.3) 



note that the exponent of N in ( |4.3|) is essentially the same as the exponent in (|3.3|) , even 
though P cnt depends on N and is one or two orders of magnitude smaller than in Section III. 



10 



FIGURES 
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Fig.l. Plot of log a max as a function of log N for the iterated quantum algorithm with minimal 
P cut = 1 - q^/(2Vn) that gtm 

improves the classical search algorithm. Even though to each N 
corresponds a different P cut , the plot still displays the universal N~ 2 / 3 dependence. 

V. CONCLUSIONS 

At the moment, quantum correction codes are restricted to the case when only a enough 
small subset of the quantum computer's g-bits undergo errors, and the probability of occur- 
rence of an error is smaller than some bound, but it is believed that quantum computations 
will be possible with physically realistic levels of noise even if the quantum correction codes 
employed undergo errors themselves. 

With this in mind, we studied the intrinsic robustness of Grover's quantum search algo- 
rithm in a noisy environment. We modelled the noise with a single parameter, the width 
of a Gaussian distribution, and allowed for two independent noises at each step of Grover's 
quantum algorithm. 
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We found that the quantum search algorithm still reaches the maximum likelihood of 
finding the searched-for state in zy/N /4 steps. The strongest effect of noise is to decrease 
the maximum probability from virtually 1 (the noiseless case) to lower values, depending on 
the size of the noise, equations ( |3.3| ) and (|3.5|). How much noise can we add to the quantum 
computer, with the criterion that a repeated application of the quantum algorithm is still 
faster than the classical one is given by equation (|4.2j ). In both cases, the allowed maximal 
noise decreases with the size of the phone book approximately as iV~ 2 / 3 . 

The presence of noise and the absence of quantum correction codes is not completely 
disastrous: the quantum search algorithm can handle by itself a reasonable amount of noise. 
Nevertheless, for large enough databases, the allowed noise becomes tiny. 

Acknowledgments. This work is supported in part by CONACYT 25504-E and 
DGAPA-UNAM IN103997. B.P.N, enjoys a scholarship from CONACYT. 

VI. APPENDIX 

The computer program we used to derive the results in section III needs an initial value of 
o (which we set to zero), and then computes P{m). If P(m max ) > P cut , then the algorithm 
increases a, repeating the process until the bound is surpassed. This gives one value for 
Cmax- We repeat the whole story again and again and average over the values of cr max found. 

Let us illustrate our procedure with an example. 

Let P cut = 0.7. For each n from 10 to 16, the program increases the value of a starting 
from in steps of da = 0.0001. The average maximal values of a thus found (in 200 runs) is 
then <7 max , shown below with its statistical uncertainty. Note that the error seems dominated 
by the step size: 
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N C^max 

1024 0.0033 ± 0.0007 
2048 0.0022 ± 0.0005 
4096 0.0014 ±0.0004 

(6.1) 

8192 0.00098 ± 0.00028 
16384 0.00070 ± 0.00020 
32768 0.00048 ± 0.00019 
65536 0.00011 ±0.00017 

Taking a smaller step, da = 0.00001, we carry through the same computations and find 
instead: 

1024 0.0022 ±0.0003 
2048 0.0015 ±0.0002 
4096 0.00095 ±0.00015 

(6.2) 

8192 0.00060 ±0.00011 
16384 0.00040 ± 0.00008 
32768 0.00026 ± 0.00006 
65536 0.00017 ±0.00005 

Curiously, when we decrease the step both the error and the central value of <r max de- 
crease. This can be understood easily, since we take as value for maximal a in each run the 
first a for which the probability after m max iterations is too small (smaller than .7 in this 
example), and thus we clearly underestimate it in gross dependence with the step. We are 
thus forced to repeat the computation of cr max and A<7 max with smaller and smaller steps, 
from d<7 = 10~ 4 to da = 10 -8 . We must now fit the dependence of <r max on da (see Fig. 2) 
and extrapolate to da = 0. 
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Fig. 2. Plot of <7 max as a function of log da: the maximum allowed value of noise characterized 
by a max before P max < -Pcut depends on the size of the step da by which a is increased in the 
program. This plot is for N = 32768 and P cut = 0.7. 

The generic relation we found is 

a max (N, da) = C (AO + f (N) da a , (6.3) 

where a = 0.30 ± 0.06 is a true constant. The values of the iV-dependent ( and £ are the 
following: 
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N 


c 




1024 


0.00104 ±0.00004 


0.0240 ± 0.0046 


2048 


0.00065 ±0.00001 


0.0163 ±0.0017 


4096 


0.00038 ±0.00002 


0.0114 ±0.0027 


8192 


0.00023 ± 0.00001 


0.0086 ± 0.0020 


16384 


0.00015 ±0.00001 


0.0079 ± 0.0025 


32768 


0.00009 ± 5 x 10- 6 


0.0068 ± 0.0020 


65536 


0.00006 ±4 x lO" 6 


0.0076 ± 0.0034 



Taking the limit dcr — > 0, we obtain the final value of u max for each N at this P cut 

N ""max 

1024 0.00104 ±0.00004 

2048 0.00065 ± 0.00001 

4096 0.00038 ± 0.00002 

8192 0.00023 ± 0.00001 ' 
16384 0.00015 ±0.00001 
32768 0.00009 ± 5 x 10~ 6 
65536 0.00006 ± 4 x 10~ 6 

The above numbers are very well fit by a straight line (in log N). 

From the data ( |6.5|) , for this value of P cut = 0.7, we find finally the relation 

<Xmax (N, Pent = 0.7) = a (P C ut) 

with a = 0.138 ± 0.012, and = -0.704 ± 0.01. 
Similarly, for other values of P cu t we found : 
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Peat OL 

0.5 0.158 ±0.011 -0.687 ±0.007 

0.6 0.146 ±0.010 -0.691 ±0.008 

(6.7) 

0.7 0.138 ±0.012 -0.704 ±0.010 

0.8 0.083 ±0.006 -0.669 ± 0.008 

0.9 0.094 ±0.001 -0.724 ±0.015 



The value quoted in the text, equation ( p.4|) , is an average of these numbers. 



To establish equation (|3.5| ), we found for each P cut a table like (|6.5| ) and then, fixing N, 



we found a good linear fit, equation Q3.5|) , with the following values of j(N) and 8(N): 



N 


7 


5 


1024 


0.0024 ±0.001 


0.0020 ±0.0001 


2048 


0.0015 ±0.00005 


0.0013 ±0.00005 


4096 


0.00092 ± 0.00003 


0.00077 ± 0.00003 


8192 


0.00057 ±0.00002 


0.00048 ± 0.00002 


16384 


0.00033 ± 0.00002 


0.00026 ± 0.00002 


13768 


0.00021 ± 0.00002 


0.00017 ±0.00002 


65536 


0.00015 ±0.00001 


0.00013 ±0.00001 



(6.8) 
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